Power-Law Sensitivity to Initial Conditions within a Logistic-like Family of Maps: 

Fractality and Nonextensivity 

U.M.S. Costa, M.L.Lyra, 
Departamento de Fisica, Universidade Federal de Alagoas, Maceio-AL, Brazil 

A.R. Plastino and C. Tsallis 

Centra Brasileiro de Pesquisas Fisicas 
Rua Xavier Sigaud 150, 22290-180 - Rio de Janeiro - RJ, Brazil 
e-mail: tsallis@cat. cbpf. br 

Power-law sensitivity to initial conditions, characterizing the behaviour of dynamical systems at 

their critical points (where the standard Liapunov exponent vanishes) , is studied in connection with 

the family of nonlinear ID logistic-like maps xt+i = 1 — a \x t \ z , [z > 1; < a < 2; t = 0, 1, 2, . . .). 

The main ingredient of our approach is the generalized deviation law limAa;(o)->o axvo) ~ H + — 

i . i 

q)\ q t\ 1 -i (equal to e 1 for q = 1, and proportional, for large t, to t 1 -i for q ^ 1; q € H is 

the entropic index appearing in the recently introduced nonextensive generalized statistics). The 

relation between the parameter q and the fractal dimension df of the onset-to-chaos attractor is 

revealed: q appears to monotonically decrease from 1 (Boltzmann-Gibbs, extensive, limit) to — oo 

when df varies from 1 (nonfractal, ergodic-like, limit) to zero. 
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I. INTRODUCTION 

The standard thermostatistical formalism of 
Boltzmann-Gibbs (BG) constitutes one of the most suc- 
cessful paradigms of theoretical physics. It provides the 
link between microscopic dynamics and the macroscopic 
properties of matter. Inspired in Shannon's Information 
Theory jjj, Jaynes reformulation of the BG theory 
greatly increased its power and scope. Jaynes provided a 
general prescription for the construction of a probability 
distribution /(x) (x £ R d stands for a point in the rel- 
evant phase space), when the only available information 
about the system are the mean values of M quantities 

(A-(x)> = J A r (x)/(x)dx, (r = l,...,M). (1) 

According to Jaynes, the least biased distribution com- 
patible with the data is the one that maximizes Shan- 
non's Information, 



Si = - J /(x) ln/(x)dx, 



(2) 



(the use of the subindex 1 will become transparent later 
on) under the constraints imposed by the mean values 
([l|) and appropriate normalization 



J /(x)dx = 1. 



(3) 



The well known answer to the above variational prob- 
lem is provided by the maximum entropy (ME) distribu- 
tion 



1 ( M \ 

foW = Y[ CXP ( X rA r (x)\ , 



(4) 



where {A r } are the M Lagrange multipliers associated 
with the known mean values, and the partition function 
Z\ is given by 



Zi 



exp 



M \ 

2J A r A r (x) dx. 

r=l / 



(5) 



Jaynes' prescription can be regarded as a mathematical 
formulation of the celebrated "Occam's Razor" principle. 
In order to obtain a statistical description of a system, 
given by the distribution /(x), we must employ all and 
only the disponible data (□), without assuming any fur- 
ther information we do not actually have. 

Jaynes informational approach allows to consider more 
general statistical ensembles than the Gibbs microcanon- 
ical, canonical, and macrocanonical. Also it provides a 
natural way to treat noncquilibrium situations. 

Despite its great success, the Boltzmann-Gibbs- Jaynes 
formalism is unable to deal with a variety of interesting 
physical problems such as the thermodynamics of self- 
gravitating systems, some anomalous diffusion phenom- 
ena, Levy flights and distributions, turbulence, among 
others (see || for a more detailed list). In order to 
deal with these difficulties, Jaynes approach is compati- 
ble with exploring the possibility of building up a ther- 
mostatistics based upon an entropy functional different 
from the usual logarithmic entropy. Recently one of us 
introduced Q the following generalized, nonextensive en- 
tropy form 



S q = 



1 



/[/(x)]*dx 



(6) 



1 



where q is a real parameter characterizing the entropy 
functional S q . This entropy recovers Si as the q = 1 
particular instance and was introduced in order to de- 
scribe systems where nonextensivity plays an important 
role; if A and B are two independent systems (in the 
sense that the probabilities associated with A + B fac- 
torize into those of A and B) we straightforwardly verify 
that S q (A + B) = S q {A) + S q {B) + (1 - q) S q (A) S q (B). 
Indeed, nonextensive behaviour is the common feature 
among the above listed problems where the usual statis- 
tics fails. The generalized nonextensive thermostatistics 
has already been applied to astrophysical self-gravitating 
systems ||, the solar neutrino problem ||, distribution 
of peculiar velocities of galaxy clusters @ , cosmology 
two-dimensional turbulence in pure-electron plasma 
anomalous diffusions of the Levy and correlated Jl I 
types, long-range magnetic and Lennard- Jones-like sys- 
tems Jl2|, simulated annealing and other optimization 
techniques Jl|, dynamical linear response theory [fl4|| , 
among others. 

The nonextensivity effects displayed by the above 
listed systems can arise from long-range interactions, 
long-range microscopic memory, or fractal space-time 
constraints. Even for dynamical systems that "live" in an 
euclidean (nonfractal) space, if the subset (of this space) 
that the system visits (most of the time) during its evo- 
lution has a fractal geometry, the generalized thermo- 
statistics might provide a better account of the situation 
than that provided by the usual statistics. Indeed, it 
is well known that nonlinear chaotic dynamical systems 
may have fractal attractors |i~5| . Two of the most impor- 
tant dynamical quantities usually employed in order to 
characterize such chaotic systems, are the Liapunov ex- 
ponents j[(|[l7), and the Kolmogorov-Sinai (KS) entropy 
p8| . In a recent effort (TPZ from here on), gener- 
alizations for these quantities inspired in the generalized 
nonextensive entropy S q (and its consequences) were in- 
troduced. The generalized Liapunov exponent X q and 
generalized KS-entropy K q provide a useful characteri- 
zation of the dynamics corresponding to critical points 
where the usual Liapunov exponent vanishes. For these 
critical cases, the exponential sensitivity to initial con- 
ditions is replaced by a power-law one, and the van- 
ishing (standard) Liapunov exponent Ai provides but a 
poor description of the concomitant dynamics. On the 
contrary, the generalized exponent X q appropriately dis- 
criminates between the different possible power-law be- 
haviours. TPZ illustrates these ideas with the (good old) 
logistic map. It is of interest to explore this formalism 
as applied to other nonlinear dynamical systems. In par- 
ticular, it is of importance to study families of dynami- 
cal systems characterized by a set of parameters. Each 
member of the family will have a different onset-to-chaos 
critical point, with a corresponding attractor character- 
ized by a Hausdorff fractal dimension df and a suitable 
value of the entropic parameter q. The study of these 
families will enlighten the relation between q and df. 

The specific aim of the present paper is to study the 



relation between the fractal dimension df of the onset-to- 
chaos attractor and the parameter q for the logistic-like 
family of maps (see [p0|-p2[ and references therein) 
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a \x t \ z , 



(J) 



(z > 1;0 < a < 2;t = 0, 1, 2, . . . ; x t e [-1,1]). Note 
that in the particular case z = 2 we recover the standard 
logistic map (in its centered representation). 

The paper is organized as follows. In Section II we 
briefly review the ^-generalizations of the Liapunov ex- 
ponent and the KS-entropy. In Section III the results for 
the logistic-like maps are presented. Our main conclu- 
sions are drawn in Section IV. 



II. GENERALIZED LIAPUNOV EXPONENT 
AND KS-ENTROPY. 

Let us consider, for a one dimensional dynamical sys- 
tem, two nearby orbits whose initial conditions differ by 
the small quantity Ax(0). We will assume that the time 
dependence of the distance between both orbits is given 
by the ansatz Jl9fl 



lim 



Ax(t) 



Ax(0)^0 Ax(0) 



[1 + (1 - q) X q t}~ (qell) (8) 



where X q is our generalized Liapunov exponent, and q is 
a real parameter characterizing the behaviour of the sys- 
tem. We verify that this equation is identically satisfied 
for t = (Vg), and that q ^ 1 yields, for large times, the 
power-law 



Ax(t) 



lim 



[{l~q)X q }— t~ (t 



(9) 



On the other hand, it is plain that for q — > 1 we recover 
the standard exponential deviation law 

Ax(t) 

lim ; ( = expfAi t] (10) 

where Ai is just the usual Liapunov exponent. The q = 1 
scenario corresponds to situations with Ai ^ 0. These 
cases describe chaotic behaviour (Ai > 0) and regu- 
lar behaviour (Ai < 0). The generalized exponent X q 
is intended to provide a convenient description of the 
marginal situations where the usual Liapunov exponent 
vanishes (Ai =0). In these last cases, we have the power 
law sensibility to initial conditions given by Eq.(^) in- 
stead of the usual exponential one. The generalized de- 
viation law (Eq.(||)) is inspired in the form of the q- 
generalized nonextensive canonical distribution, given by 



Pi = 



[1 - (l-gJ/JcjVd-d 



z„ 



(11) 



2 



with the generalized partition function being given by 



that under the assumption of equiprobability reduces to 



1/(1-9) 



(12) 



where (3 = 1/kT and {tj} is the full set of eigenvalues of 
the Hamiltonian of the system. Notice that, in the limit 
q — * 1, this thermal canonical equilibrium distribution 
reduces to the ordinary BG one 



Pi 



exp[-/3ej] 



with 



Zi = ^cxp[-/3e;] 



(13) 



(14) 



It is worth to remark that the marginal case with van- 
ishing (standard) Liapunov exponent Ai displays a very 
rich and complex behaviour, reminiscent of what happens 
at the critical point of thermal equilibrium critical phe- 
nomena. To just say that Ai = is a very poor descrip- 
tion of its richness, intimately connected to fractality. 
Indeed, within our generalized formalism, the parameter 
q provides a characterization of the kind of power-law 
sensitivity to initial conditions involved, and is expected 
to be related to the fractal dimension df of the corre- 
sponding attractor. 

In order to discuss the generalized KS-entropy, let us 
consider a partition of phase space in cells with size char- 
acterized by a linear scale I. We will study the evolution 
of an ensemble of identical copies of our system. We as- 
sume that all the members of the ensemble start at t = 
with initial conditions belonging to one and the same 
cell. This means that the probability associated to that 
privileged cell is 1, while the remaining cells of the par- 
tition have vanishing initial probabilities. As time goes 
by, and due to the sensitivity to initial conditions, our 
ensemble will spread over an increasing number of cells. 
The standard KS-entropy can be regarded as the rate of 
growth of the Boltzmann-Gibbs entropy associated with 
the partition probability distribution. 

Within the generalized nonextensive thermostatistics, 
the entropy functional for a discrete probability distribu- 
tion {pi} is given by 



1 



which, for equiprobability, becomes 

w x -i - 1 



1-9 



(15) 



(16) 



The use of Eq. (|l5|), instead of Si = — J^ i=1 p% fripi, 
yields (along Zanette's lines p3|) to the following gener- 
alization of the Kolmogorov- Sinai entropy 



K, 



q = lim km Jim — (S q (N) - 5,(0)) , 

t— »0 I— >0 N— >oo JVT 



(17) 



A" 



lim 



lim lim _. 

t->0£->0 N^oo Nt 



1 [W{N)\ 



(1-9) 



1 - 



(18) 



In both equations (|IJJT^) we have maintained the tradi- 
tional t — > which applies for a continuous time t\ it is 
clear however that, in our present case, this limit does 
not apply since our t is discrete. We must remark that 
our generalization K q of the KS entropy is different from 
the generalizations A(/3) based upon Renyi informations, 
usually called " Renyi entropies" in the literature of ther- 
modynamics of chaotic systems [ p4[ (sometimes the pa- 
rameter characterizing these generalizations of the KS 
entropies is called q instead of (3 p5[ . This parameter q 
should not be confused with our q). 

Consistently with the behavior indicated in Eq. (^), we 
have (along Hilborn's lines |17j]) 



W{N)= [! + (!-<?) \ q Nt]- 



which, replaced into Eq. 
dynamical systems) 



(19) 



8[), immediately yields (for ID 



K„ = A„ 



(20) 



This relation holds if A, > (K q vanishes if A, < 0); 
it constitutes a generalization of the well known Pesin 
equality Ai = Ai (if Ai > 0; Ai = otherwise), and 
unifies (within a single scenario for both exponential and 
power-law sensitivities to initial conditions) the connec- 
tion between sensitivity and rythm of loss of information. 



III. THE LOGISTIC-LIKE MAP. 



Let us now illustrate some of the above concepts by 
focusing the logistic- like maps ([?]) . These maps are rela- 
tively well known and have been addressed in various oc- 
casions ( pl| , p2| and references therein). The topological 
properties associated with them (such as the sequence of 
attractors while varying the parameter a) do not depend 
on z, but the metrical properties (such as Feigenbaum's 
exponents) do depend on z. We shall exhibit herein that 
the same occurs with q. Indeed, although quite a lot is 
known for these maps, their sensitivity to the initial con- 
ditions at the onset-to-chaos has never been addressed 
as far as we know. As we shall see, for all values of z, 
the sensitivity is of the weak type p9[ , i.e., power-laws 
instead of the usual exponential ones. 
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FIG. 1. Log-log plot of lim A;c(0) ^o (xfpy) versus tne 
number of iterations N calculated for xq = (the slope 
1/(1 — g) is calculated using the upper bound points): (a) 
2 = 1.25; (b) ^ = 2; (c) z = 3 



We present now our main numerical results. We com- 
puted, as functions of z, the parameter q and the critical 
fractal dimension df (determined within the box count- 
ing procedure). 

In Fig. 1 we exhibit, for typical values of z at its 
chaotic threshold a c (z) and using xo — 0, a plot of 
lnlim Ax(0 )^Q = J2^=i iia [ az \ x n\ Z ~ 1 } versus In N, 

where N is the number of iterations. (Notice that for 
convenience we use, as argument of the logarithm, not 
exactly the derivative dxt+i/dxt , but rather its absolute 
value). For each of these plots we see an upper bound 
whose slope equals 1/(1 — q) (see Eq. (^)), from which 
we determine q. 
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FIG. 2. z-dependence of the entropic index q. The inset 
conveniently represents the z — > oo neiborhood (the continu- 
ous line is the best fitting with a curve q = 1 — ao/(z — l) ai ; 
we obtained ao — 0.746 and a\ — 0.613). 



In Fig. 2 we show the behaviour of the parameter q 
as a function of the parameter z characterizing the map. 
The figure suggests that for z — > 1, q tends to — oo, while 
in the limit z — > oo, q approaches 1. In Fig. 3 we can 
see, for typical values of z at its chaotic threshold a c (z), 
the number of filled boxes as function of the number of 
boxes, corresponding to the box counting method em- 
ployed in order to determine the fractal dimension df. In 
Fig. 4 the behaviour of the fractal dimension df of the 
chaotic critical attractor as function of z is depicted. We 
can observe that, as z — * 1, df seems to go to 0, while, 
in the limit z — > oo, the c?/-curve approaches unity.In 
Fig. 5 we show the behaviour of the parameter q as a 
function of the fractal dimension df. We can see that q 
displays a monotonically increasing behaviour with the 
fractal dimension df. 
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FIG. 3. Box counting procedure for determining the on- 
set-to-chaos fractal dimension df for typical values of z. 
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FIG. 4. z-dependence of the fractal dimension df. 
The inset conveniently represents the z — * oo neibor- 
hood (the continuous line is the best fitting with a curve 



bi ■■ 



= exp - 
0.27). 



[bo/(z — 1) 1 ]; we obtained bo 



0.62 and 



z 




(7 




1 


1* 


— OO* 


0* 


1.05 


1.0816488... 


-4.52 ± 0.03 


0.24 ±0.02 


1.10 


1.1249885... 


-2.28 ±0.02 


0.32 ±0.02 


1.25 


1.2095137... 


-0.76 ± 0.01 


0.40 ±0.01 


1.5 


1.2955099... 


-0.12 ± 0.01 


0.48 ± 0.01 


1.609 


1.3236435... 


0.00 ± 0.01 


0.50 ±0.01 


1.75 


1.3550607... 


0.13 ± 0.01 


0.52 ± 0.01 


2.0 


1.4011551... 


0.24 ±0.01 


0.54 ±0.01 


2.5 


1.4705500... 


0.41 ±0.01 


0.57 ±0.01 


3.0 


1.5218787... 


0.47 ±0.01 


0.60 ±0.01 


5.0 


1.6456203... 


0.69 ±0.01 


0.65 ±0.01 


oo 


2* 


1* 


1* 



TABLE I - * denotes the limiting values suggested by the 
numerical results; a c = 1 for z = 1 also has analytic sup- 
port (see pl| and references therein); a c — 2 for z — > oo 
also has renormalization group support pM. 



IV. CONCLUSIONS. 



The particular value q — 0, that describes linear sen- 
sitivity to initial conditions corresponds, with notable 
numerical accuracy, to the fractal dimension df = 0.5 
(numerically df = 0.50 ± 0.01, occuring for z = 1.609 ± 
0.001). It is remarkable that as the fractal dimension 
tends towards 1, the parameter q approaches 1. If this 
tendency becomes confirmed by analytic results or more 
powerful numerical work, this would be very enlighten- 
ing, because in the limit df — ► 1 the attractor loses its 
fractal nature (in the sense that df coincides with the 
euclidean dimension d = 1), and the usual statistics (i.e., 
the usual exponential deviation of nearby trajectories), 
characterized by q — 1, would be recovered. On the 
other extreme, as df — > 0, q appears to approach — oo, 
hence 1/(1 — q) = 0, which can be considered as an indi- 
cation of a possible logarithmic sensitivity to the initial 
conditions. Our results are summarized in Table I. 
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FIG. 5. d/-dependence of the entropic index q; the inset 
conveniently represents the df < 1/2 region. Notice that the 
point (df,q) — (1/2,0) seems to belong to the curve. The 
dashed lines are guides to the eye. 



We have exhibited, for a family of logistic-like maps 
, the behaviours of the entropic parameter q and the 
fractal dimension df of the onset-to-chaos attractor. We 
showed that, at this critical point, power deviation laws 
for nearby orbits, similar to the ones appearing |fl9|| in 
the logistic map, are observed. The concomitant value 
of q is related to the chaotic attractor fractal dimen- 
sion. It would no doubt be interesting to find out if, 
for generic nonlinear dynamical systems, q depends only 
on df (support of the visiting frequency function) or 
also upon other characteristics of the critical attractor, 
such as the visiting frequency function itself. In order 
to answer this question, it would be useful to explore 
the behavior of families of maps whose possible chaotic 
critical points depend on more than one parameter. Such 
studies, as well as the application of these concepts to 
self-organized criticality [26j, would be very welcome. 
Finally, let us stress that the present study provides a 
direct and important insight onto a problem which has 
been quite cllusive up to now, namely the microscopic in- 
terpretation of the entropic index q characterizing nonex- 
tensive statistics. The present results clearly exhibit that 
what determines q is not the entire phase space within 
which the system is allowed to evolve (the euclidean in- 
terval — 1 < Xt < 1 in the present examples), but the 
(possibly fractal) subset of it onto which the system is 
driven by its own dynamics. Consistently, whenever the 
relevant fractal dimension approaches its associated eu- 
clidean value (d = 1 in the present case), extensivity 
(i.e., q = 1) and standard BG thermostatistics naturally 
become, as is well known, the appropriate standpoints. 
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